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We introduce a new version of non-linear electrodynamics which is produced by a spontaneous 
symmetry breaking of Lorentz invariance induced by the non-zero expectation value of the gauge 
invariant electromagnetic field strength. The symmetry breaking potential is argued to effectively 
arise from the integration of massive gauge bosons and fermions in an underlying fundamental 
theory. All possible choices of the vacuum lead only to the remaining invariant subgroups T(2) 
and HOM(2). We explore in detail the plane wave solutions of the linearized sector of the model 
for an arbitrary vacuum. They present two types of dispersion relations. One corresponds to the 
case of the usual Maxwell electrodynamics with the standard polarization properties of the fields. 
The other dispersion relation involves anisotropies determined by the structure of the vacuum. The 
corresponding fields reflect these anisotropies. The model is stable in the small Lorentz invariance 
violation approximation. We have also embedded our model in the photon sector of the Standard 
Model Extension, in order to translate the many bounds obtained in the latter into corresponding 
f"*^ ' limits for our parameters. The one-way anisotropic speed of light is calculated for a general vacuum 

and its isotropic component is strongly bounded by 8c/ c < 2 x 10~ 32 . The anisotropic violation 
contribution is estimated by introducing an alternative definition for the difference of the two- 
way speed of light in perpendicular directions, Ac, that is relevant to Michelson-Morley type of 
experiments and which turns out to be also strongly bounded by Ac/c < 10 -32 . Finally, we 
. speculate on the relation of the vacuum energy of the model with the cosmological constant and 

O l' propose a connection between the vacuum fields and the intergalactic magnetic fields. 
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I. INTRODUCTION 



The possible violation of Lorentz invariance has recently received a lot of attention, both from the experimental 
\ and theoretical sides. In the latter case, mainly in connection with possible effects arising from drastic modifications 
of space-time at distances of the order of Planck length suggested by most of the current quantum gravity approaches. 
Experiments and astrophysical observations set stringent bounds upon the parameters describing such violations, 
which nevertheless are still been subject to improvement in their precision. Broadly speaking, there are two basic 
t-H ' possibilities which produce such a breaking: (1) one introduces by hand a number of non-dynamical tensor fields, 
whose fixed directions induce the corresponding breaking in a given reference frame. Examples of this approach 
are the phenomenological Myers-Pospelov model [l| together with QED in a constant axial vector background @|. 
(2) a second possibility is to dwell upon spontaneous symmetry breaking (SSB). This is basically what is done in 
the Standard Model Extension Q, [j], where such non-dynamical tensor fields are assumed to arise from vacuum 
expectation values of some basic fields belonging to a more fundamental model, like string theory for example. 

Nevertheless, we are interested in studying a model of electrodynamics which incorporates spontaneous Lorentz 
symmetry breaking, without having to go to a string theory setting. In order to produce such a model we require the 
presence of the non-zero vacuum expectation value of a tensorial field of rank greater than zero. 

A first possibility arises from having a non-zero VEV of the potential field A^. This case has been thoroughly 
studied and leads to the interesting idea that the photon arises as the corresponding Goldstone boson of the global 
spontaneous Lorentz symmetry breaking. Normally the masslessness of the photon is explained by invoking abelian 
gauge invariance under the group U(l). This almost sacred principle has undoubtedly been fundamental in the 
development of physics and it is very interesting to explore the possibility that it could have a dynamical origin 
Q . This idea goes back to the works of Nambu and Biorken Q together with many other contributions [1, @] . 
Recently it has been revived in references [l(J [HI 0, EH 13 ■ One of the most explored approaches along these lines 
starts form a theory with a vector field endowed with the standard electrodynamics kinetic term plus a potential 
designed to break Lorentz symmetry via a non-zero vacuum expectation value < >, which defines a preferred 
direction in space-time. This potential also breaks gauge invariance. These are the so called bumblebee models [15| . 
The subsequent symmetry breaking, obtained from the non-zero minimum of the potential, splits the original four 
degrees of freedom in _B M into three vectorial Nambu-Goldstone bosons A^, satisfying the constraint A^A^ — ±M 2 , 
to be identified with the photon, plus a massive scalar field a, which is assumed to be excited at very high energies. 
Under this threshold one basically recovers the Lagrangian for electrodynamics in vacuum plus the above mentioned 
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non-linear constraint, which is interpreted as a gauge fixing condition. Calculations at the tree level [6| and to the 
one loop level [l2j have been carried on, showing that the possible Lorentz violating effects are not present in the 
physical observables and that the results are completely consistent with the standard gauge invariant electrodynamics. 
This is certainly a surprising result, which can be really appreciated from the complexity of the calculation in Ref. 

Nevertheless, alternative kinetic terms in the bumblebee models can lead to a theory differing substantially from 
electromagnetism [ltij . The use of nonpolynomial interactions within this framework has been explored in [17] . The 
idea of a photon as a Goldstone boson has also been extended to gravitons in general relativity [T^, [TH, [l{| . 

It is noteworthy to recall here that the constraint A^A^ = ±M 2 was originally proposed by Dirac as a way to 
derive the electromagnetic current from the additional excitations of the photon field, which now ceased to be gauge 
de gree s of freedom, avoiding in this way the problems arising from considering point-like charges [20]. As reported 
in [21] , this theory requires the existence of an ether emerging from quantum fluctuations, but not necessarily of the 
violation of Lorentz invariance, due to an averaging process of the quantum states producing such fluctuations. 

In this work we explore the possibility of having a non-zero VEV of the electromagnetic strength, so that we avoid 
the problem of disturbing gauge invariance from the very beginning. To this end we consider a non-linear version of 
Electrodynamics with a potential term having an stable minimum for a constant electromagnetic strength. The paper 
is organized as follows: In Section II we make it plausible the existence of such gauge invariant potential V e f f(F 2 ), 
arising as an effective low energy contribution from the integration of degrees of freedom corresponding to massive 
gauge-bosons and fermions in a more fundamental theory. This potential also exhibits the right behavior in the low- 
intensity field approximation. Section III contains the construction of the model starting from a non-linear version of 
Electrodynamics which is spontaneously broken by choosing a vacuum characterized by a constant electromagnetic 
tensor . In Section IV we discuss the possible realizations of such vacuum in terms of the associated electromagnetic 
fields e and b, identifying the Lorentz subgroups that remain invariant after the breaking. The equations of motion 
corresponding to the propagating sector, i.e. those arising from the quadratic terms in the Lagrangian, are examined 
in Section V using a covariant approach, as well as the standard 3 + 1 decomposition. Also, the modified Maxwell 
equations are presented. Section VI contains a summary of the modified photon dispersion relations u> — w(k), 
together with the expressions for the electromagnetic fields E and B of plane waves propagating in the different vacua 
parameterized by e and b. The model is shown to be stable in the small Lorentz invariance violation approximation. 
In section VII we study the embedding of our model in the corresponding sector of the Standard Model Extension 
[23 | and use the experimentally established bounds upon the LIV parameters to constrain the electromagnetic fields e 
and b. Finally we conclude with a summary of the work together with some comments in Section VIII. The Appendix 
contains some details of the calculation leading to the motivation described in Section II. 



II. MOTIVATION 



One of the main challenges of the proposed approach is to make it plausible the existence of a gauge invariant 
potential V(F a p) depending upon the electromagnetic tensor, which is stable for large values of the electric and 
magnetic fields and which possess a minimum for some constant value of F a p. In particular, we will consider at this 
level only the case of magnetic fields in order to avoid further instabilities due to pair creation in strong electric fields. 

A first possibility that arises is to consider the effective photon interaction in QED after integrating the fermions. 

Unfortunately, a detailed calculation of this potential, in the one-loop ( Vg^f ) and two-loop (v^j^j approximations 
[25l | . shows that it is given by 

V EM {B) = \b 2 + + V<& = \b 2 - ^- 2 B 2 

in the limit when B >> Bq = m 2 /e, which clearly is unbounded from below. 

A new possibility arises in a further generalization of the mechanism proposed in Refs. 0,0]! restricted to the gauge 
invariant case. We start from a conventional gauge theory including fermions, gauge fields and Higgs fields which 
provide masses, via spontaneous symmetry breaking, to the gauge fields, except for the photon potential A^, which 
carries the electromagnetic U(l) gauge invariance. For simplicity we consider only one fermion species, as in Ref.Q, 
but the idea can be generalized to many of them. In the latter case, the subsequent integration over the fermionic 
degrees of freedom could be made only for those having a mass larger than a given scale, thus leaving unintegrated 
the lightest ones. 

Integrating the massive gauge bosons of this theory leads to a power series in bilinears in the fermion fields, which 
are suppressed by powers of the energy scale A that characterizes the SSB mechanism providing the mass to the gauge 
vector bosons [8| . We focus upon the neutral sector of the model where the following contributions to the effective 
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We are considering all possible contributions which are consistent with gauge invariance, the remaining symmetry 
of this sector of the theory. Notice also that terms containing the U(l) covariant derivative D v = [d u + ieA v ) are 
further suppressed with respect to those without it. Also we consider only the n = 1 case, which provides a gauge 
invariant generalization of the model in Ref. 7]. Thus our starting point is 
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The index M : S,V,T, PV, PS labels the tensorial objects that constitute the Dirac basis: Scalar, Vector, Tensor, 
Pseudovector and Pseudoscalar, respectively. The generic index a labels the covariant (contravariant) components 
of each tensor class. More details are given in Table I of the Appendix. We follow the Dirac algebra conventions of 
Ref. dU, with appropriate factors chosen in such a way that each current (^M ^) is real. In this way, the proposed 
generalization includes all possible fermionic quartic interactions. Also, the coefficients tm are assumed to be of order 
one. 

The fermions can be integrated next by introducing the corresponding auxiliary fields (Cm) for each real current 
(jm)q = ^M a ^, as suggested in Ref. [8J, and following standard manipulations in the path integral formulation. 
Some details are given in the Appendix. 

Our main goal is to estimate the additional corrections provided by such currents to the purely electromagnetic 
effective potential (|T|) and to investigate whether or not it is possible to generate a net positive contribution in the 
high- intensity field limit. To have a preliminary simpler expression for the resulting extension we choose to consider 
separately each of the contributions arising from the different currents and only add them to the total effective potential 
at the end. In this way we are nor considering possible interference terms among the currents, which nevertheless 
could be calculated within the approximations performed. Calling V^ 2 ' the finite part of the total effective potential 
arising from the currents j'mj given in Eq. (|259p of the Appendix, we obtain 
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Here m is the mass of the integrated fermion which is such that m > A. The main point of the estimation is that we 
obtain 



= tm 0m = 4r^ + 16rv + 48r^ — 16rpv — 4rpg, 



(6) 
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which can be made positive from a judicious choice of the parameters tm- 

In the high-intensity field approximation, the dominant term in the electromagnetic contribution goes like 
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which can be overcame by the contribution from the additional currents 
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provided 6 > 0. 

As we mention in the Appendix, this preliminary estimation has only considered an appropriate regularization 
of the divergent integrals, but is missing an adequate renormalization of the total effective Lagrangian L e ff. Our 
renormalization conditions would be 



lim L 



eff 



P B 2 



p>0, 



(9) 



4 



in such a way that the effective potential is a decreasing function in the vicinity of B = 0. The final normalization to 
the value of £ e // = —b 2 /2 will be imposed at the end of the procedure, after we have expanded the magnetic field 
around an stable minimum B^in of V e t / , such that 

B = B Mm + b (10) 

and we are left with the physical component b. The renormalization procedure should be analogous to the one carried 
in Ref. [25| for the two- loop contribution L^ M (A) to the effective electromagnetic action and it is rather involved due 
to the complicated field dependence of the regularized version of the divergent integrals. This calculation is beyond 
the scope of these preliminary estimations. 

Having an effective potential with the characteristics described above guarantees the existence of an absolute 
minimum, which we assume to arise at values lower that the critical magnetic and electric fields, in such a way to 
avoid fermion pair production in the later case. We also generalize this idea to the case of an arbitrary constant 
electromagnetic field F^ v . The study of the stability of the model under radiative corrections is beyond the scope of 
the present work. 



III. THE MODEL 



In the previous section we have made it plausible the existence of a gauge invariant potential V e ff(F 2 ) having an 
stable minimum and arising as an effective low energy model from the integration of degrees of freedom corresponding 
to massive gauge-bosons and fermions in a more fundamental theory. This potential also exhibits the right behavior 
in the low-intensity field approximation. To examine the dynamical consequences of the symmetry breaking we have 
to adopt a convenient parametrization of such potential. There are various choices in the literature which would 
translate into the following for our case [23J 

V(F) = X(F 2 ± C 2 ), V(F) = -^{F 2 ± C 2 ) 2 , (11) 

where is the constant value of the electromagnetic field at the minimum. Instead we will adopt the standard 
Ginzburg-Landau parametrization 

V{F^)^ l -aF 2 ^^{F 2 ) 2 . (12) 

Since the basic variable we want to deal with is F^, we start from the effective Lagrangian 

L(F a0 , X„) = -V{F af} ) - F^d v X^ F af} = -F Pa , (13) 

where the dual field F v ^ is 

F Ufi = ^e v ^F aP (14) 

and the are just Lagrange multipliers that will finally impose the condition that the field strength satis- 
fies the Bianchi identity and can then be expressed in terms of the vector potential. The above Lagrangian 
is similar to that used in the standard formulation of non-linear electrodynamics [26| . Our conventions are 
= diag{+, -,-,-), e 0123 = +1, e 123 = +1. 
The next step is to determine the vacuum configuration of the theory, corresponding to the minimum of the energy 
E. The energy density T 00 is calculated via the standard Noether theorem starting from the Lagrangian (Ti"3"|) and 
produces 



E = / d A x T w = / d 6 x [V(F al3 ) + F la (da*)] ■ (15) 



The extremum conditions, obtained by varying the independent fields F a p and X a , are 

SE dV 
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In order to maintain four dimensional translational invariance we look for extrema of the fields which are independent 
of the coordinates. From Eqs. p^|) and (JT7J) we verify that this is possible provided the condition 



dV 

dF pa j Extr 



= 0. (18) 



holds. In other words, to find the vacuum configuration we have to extremize the effective action, subjected to the 
condition that F a p and X a are constant fields. Applying these requirements to (TT5|) plus the choice (|T2l we obtain 

dV 

w - = 0=( a + ( 3F")F^ (19) 
which is solved by a constant (FExtr) a p = C a p such that 



(f\ = ~ = C^0. (20) 

The expansion around the minimum (C^, Cp) is written as 

F a0 (x) = C af3 + aap(x), X ll = Cp+ X„. (21) 
Introducing such expansion in the potential (|12p , we arrive to 

V(F al3 ) = V{C a!3 )+(3{C a!3 a^) 2 + (3(C aP a a f s )(a^a^) + ^(a a0 a a ' 3 ) 2 , (22) 
= V{C a p) + V(a aP ). (23) 

We notice that (|22|) has no linear term in a ai , and also that the quadratic term j3 (C Q( 3a Q7 ) 2 is positive provided 

f3 > 0, (24) 

thus verifying the expansion around a minimum. On the other hand, according to Eq. (|20p. the sign of a is arbitrary, 
in such a way that it is opposite to the sign of C 2 . 

In this way, the spontaneously symmetry broken action is 



S(a a p, Xp) = - Jd 4 x (^e^ af} a aP d v Xp + V{a a p) + V(C)^j 



(25) 



where we have eliminated the total derivatives arising from the shifts of the fields J21|) performed in the original 
Lagrangian (TIB)) . 

Next we consider the equations of motion derived from (I25|) and show that the Lagrange multiplier X^ is fully 
determined up to gauge transformation X^ — > X^ + d^x- The equations are 

5a a p : - € "i^d v X lt -2^- = 0, (26) 

da a/3 

5X M : e v ^d v a a0 =Q. (27) 

Eq. (|2T|) establishes that the two form a is closed, i.e. da = 0. The Hodge-De Rham theorem says that the most 
general solution to this is obtained by requiring 

a = dA + ls, (28) 

where A is a one form, / is a constant and s is an harmonic two-form. Since in our case we naturally have one such 
form at our disposal, arising precisely from the chosen vacuum, we take 



s = -C^dxPdaf, (29) 



which is certainly harmonic because it is constant. Calling dA — f we have 

a a p(x) = lC a p + f a p{x). (30) 
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From (l26l) we obtain 

B. < 

which are the generalization of Maxwell equations for the potential A a . From (f2l>)) we obtain d„X p — d p X v yielding 



dV 

(d p X a - d a X p ) = € paaP - = R p<J (Ap). (32) 

oa a p 

We verify that R pa (Ap) can in fact be written as the left hand side of Eq. (|32|) by calculating € a ^ f " J dpR pcr . The result 
is zero in virtue of the equations of motion for Ap. Fixing the corresponding gauge freedom by setting d^X^ = 0, for 
example, we obtain 

d 2 X a = d p R pa (A ), (33) 

which completely determines the Lagrange multiplier X a . Going back to the action (125|) . integrating by parts and 
substituting the constraint (j27|) we are left with the reduced effective action 



S(A a ) = - J d A x (V(a a p) + V(C)) , a a p(x) = lC a0 + d a A p - d p A a . (34) 
The resulting equations of motion are 

5A a : dp (^-) = 0, (35) 



dap 

which coincide with the previous ones in Eq. (|31[) . 

It is convenient to rename the dimensionless parameter / = £ — 1. Substituting (|30[) in (|34[) and defining 

C„ v = ^D^, B =^>°> ( 36 ) 

where 

e= 2(e-lUC<*eC a py (€ a -l)(^^)>Q, (37) 
we obtain our final effective action 

S(A a ) = J d 4 x (- \ [1 - D 2 B] D 2 - \f p „r - B [{D^D + {f^D?) , (38) 

where f a p(x) = d a Ap — dpA a . We recognize the standard Maxwell kinetic term in the right hand side of Eq. (f55| . 
The only restriction now is B > 0, with D 2 arbitrary. 

Let us emphasize that the case of purely spontaneously broken Lorentz invariance really corresponds to the singular 
choice £ = 1,(1 = 0). In this case, the corresponding action would not possess the standard kinetic term —\j pv ]^ v but 

will start with the quadratic term (D pv j ^ v ) 2 . This very interesting case is not considered here and its investigation 
is postponed for future work. 

IV. SYMMETRY ALGEBRAS ARISING FROM DIFFERENT CHOICES OF THE VACUUM 

In this section we study the possible vacua allowed by the tensor symmetry breaking and we also identify the 
corresponding subgroups of the Lorentz group which are left invariant after the breaking. In order to do this, it is 
convenient to parameterize the background field D pv in terms of three dimensional components 

Dij = —eij m b m , D 0i = ej, £123 = I (39) 



ei e 2 e 3 

-ei —63 6 2 

-e 2 63 -bi 

-e 3 —b 2 bi 



(40) 
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which will mix when going to another reference frame via a passive Lorentz transformation. Since we have two vectors 
that determine a plane we choose a coordinate frame where 

b=(0,0,&), e=(0, e 2 = \e\smip, e 3 = | e l cos VO, (41) 

That is to say, we have chosen the plane of the two vectors as the (z — y) plane, with the vector b defining the 
z-direction and ip being the angle between b and e. In this way the matrix representing the vacuum is 






-e 2 



e 2 
-b 
b 



-e 3 



e3 






(42) 



The most general infinitesimal generator G, including Lorentz transformations plus dilatations is 



G = [G M „] = -i 



Z Xl x 2 x 3 

%i z -2/3 2/2 

X2 2/3 z -j/i 

X3 -2/2 2/1 z 



(43) 



Motivated by the work in Ref. [27j we are including conformal dilatations D among our generators. Within this 
restricted Poincare algebra, this generator commutes with the remaining ones corresponding to pure Lorentz trans- 
formations and can be realized as a multiple of the identity. We do this in order to explore the possibility of having 
the largest possible invariant sub-algebra after the symmetry breaking. 

The condition for the vacuum to be invariant under the infinitesimal transformations is 



= G M n D c 



and the resulting equations are: 



x 2 b-y 2 e 3 = -y 3 e 2l x 2 e 3 + y 2 b = x 3 e 2 , 
x\e 2 + 2zb = 0, yie 2 + 2ze 3 = 0, 
-xib + yxe 3 = 2ze 2 , xie 3 + yib = 0. 



(44) 



(45) 
(46) 
(47) 



Next we consider the solutions of the above system of equations corresponding to the seven non-trivial possibilities 
dictated by the choices in the arrangement (6, e 2 , e 3 ) in the above coordinate system. In the sequel we use the 
notation of Ref. [28] for the Lorentz group generators. 



A. Case b — e 3 = 0, e 2 / 

This leads to 

z = 0, xi = yi = 0, x 3 = y 3 = 0. 
We have two free parameters, x 2 , y 2 and the generator is 

G = x 2 K 2 - y 2 J 2 . 



(48) 
(49) 



B. Case b = 0, e 3 / 0, e 2 = 

Here we have 

z = 0, xi = yi = 0, x 2 = y 2 = 0, 
so that we are left with a two parameter (x 3l y 3 ) subalgebra with generator 

G = x 3 K 3 -y 3 J 3 . 



(50) 
(51) 



C. Case b 0, e 3 = 0, e 2 = 



This case is analogous to the previous one. We have 

z = 0, x\ = yi = 0, x 2 = 2/2 = 0, 
with a two parameter (X3, 2/3) subalgebra. 

D. Case 6 = 0, e 3 / 0, e 2 / 

In this case we have 

z(el + el) = -> 2 = 0, 

£1 = 0, yi = 0, 2/2 = 2/3 — , £2 =£3 — , 

e3 e 3 

yielding a two parameter (£3,1/3) subalgebra with generator 



Calling 



we can show that 



Here we have 

and the consistency condition 
which leads to two possibilities 

Here we have 



a= (|# 2 • a 1 ) //,(./•'■ • |j 2 ) 



G X3 = ^K 2 +K\ G y3 = - (J 3 + ^J 2 
e 3 m V e 3 



[G33, Gyg] — 0. 

E. Case 6/0, e 2 / 0, e 3 = 



2/3 = -X2-, x 3 = +2/2— ,2/1=0 
&2 e 2 



- e 2 .) = 0, 



1. Subcase z = 



6 6 

2/3 = -x 2 — , x 3 = +2/2— , 2/1 = 0, xi = 
e2 e2 



and we are left with a two-parameter subgroup where 



Calling 



G = .r 2 [K 2 + ^A+y2(^K 3 -J 2 ). 



G J2 = (K* + ±A, Gy2 = (±K*-J*). 



we can verify that 

[G X2 , G V2 ] = 0. 
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2. Subcase z / 

Here we have 

b 2 - e\ = -► b = se 2 , s = ±1, (64) 

which corresponds to a plane wave VEV. 
Summarizing we have 

2/1 = 0, x 1 = -2sz, x 2 = -sy 3 , y 2 = sx 3 , (65) 
Here we are left with a three parameter Lie algebra (z, x 3 , y 3 ) and the generator is 

G = —z (2SK 1 + il) + x 3 (K 3 ~ sJ 2 ) - y 3 (sK 2 + J 3 ) . (66) 

Defining 

G z = - (2-sK 1 + il) , G x = K 3 ~ sJ 2 , G y = - (J 3 + sK 2 ) , (67) 

we obtain the algebra 

[G Z ,G X ] = 2iG x , [G z ,G y ] = 2iG y , [G X) G v ] = 0, (68) 

which is isomorphic to HOM(2). 

The condition (|64| implies that D 2 = D^D^ — 0. Since B defined in Eq. (|36|) has to be finite and positive we 
have to be careful in the limiting procedure leading to this quantity. Also we need to have C ilv G^ v =^ 0, in order that 
the original parameters a and (3 are well defined. Let us consider 



such that 



The choice 



with 7 arbitrary, guarantees that 



is finite. Besides we must demand that 



Next we consider 



e 2 = sb + e, e 3 = 0, e , (69) 
D 2 = 2{b 2 - ej) = -Asbe. (70) 
£ = 7\/i, (71) 



6 =£ b ^ 



sb > 0. (73) 



which ensures that C 2 is finite. The fact that C 2 is negative is consistent with the second condition in (|37|) . Summa- 
rizing, in this case we have 

/3 = ^>0, a = \, (75) 
showing that the limiting procedure is well defined. 



F. Case 6/0, e 3 / 0, e 2 = 



We have 

z = 0, Xi = yi = 0, x 2 = 2/2 = 0. 

The free parameters here are x 3 ,y 3 . 



(76) 
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G. Case 6^0, e 2 # 0, e 3 / 

Here 



2/3 = -^2 1-2/2 — , £3=2:2 h 2/2 — , (77) 

e2 e 2 e 2 e 2 



together with the consistency condition 



z(b 2 -e\- ej) = 0, zbe 3 = 0, (78) 



-"3 °2 

which leads to z — 0. We have two free parameters: £2,2/2 ■ The generators are 

^2 , ^3 , 1 rsA _ / 7 2 _ A ^3 , j3 



Dehning 



G = x 2 \K 2 + ^K* + - J 3 - y 2 J" - -X 3 + ^ J a . (79) 
e 2 e 2 / V e 2 e 2 



G' X2 = (K 2 + ^K 3 + -jA, G '=-(j 2 --K 3 + ^J s ), (80) 
"- e 2 e 2 / a2 V e 2 e 2 



we can verify that 

[G' X2 , C; 2 ]=0. (81) 

Summarizing, all the two-parameter subalgebras that leave the vacuum invariant are isomorphic to T(2), while the 
only three-parameter subalgebra, corresponding to the case (E-2), is isomorphic to HOM{2). 

H. THE SYMMETRIES OF THE ACTION 

To conclude this Section, let us consider the full effective action 

l = - \f^r - b [(DVf MV ) + (u.rnf m 

and verify the invariance under the Lorentz plus scale transformations that leave the vacuum invariant given by the 
generator (|4"3")l together with the condition (|44"1) . The term fauf^" is invariant under the full Lorentz group plus scale 
transformations, which includes also the second term in the parenthesis proportional to B. Let us verify the invariance 
of the term {D^ v /^„) under the transformation (|44|) . We have 

5 (D^U U ) = (SD^) + (D^Sf^) = (D^Sf^) , (83) 

because (144)) is the subset of transformations that leave the vacuum D pv invariant. The remaining transformations 
are given by those of a contravariant tensor 

(D^Sfnu) = -D^ (fauG^ + f m G a v ) = -f p „ (G p a D aa + G a a D pa ) = (84) 

because the last equation in the right hand side just reproduces the vacuum invariance condition (|44[) . In these way, 
the symmetries of (|82|) are in fact those of the vacuum D pi/ . 



V. THE EQUATIONS OF MOTION: THE PROPAGATING SECTOR 

To study the propagation properties we consider only the quadratic terms in the effective Lagrangian 

io = -^r-BM"") 2 , (85) 
where we recall that f a p = d a Ap — dpA a . The equations are 

(d 2 A p - dpd a A a ) = -%BD af) d a (D^dM . (86) 
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We have verified the consistency of the above when taking . 
It is convenient to dehne 



X = D^dpA,, (87) 



and to introduce the notation 

In this way we have 
together with 



D ak d a = D k = [D ok d - Dikdt] , 

D = D i0 di. (89) 



X = - {DjAj + D Ao) , (90) 
D d = -diDi. (91) 



A. Covariant formulation 

The modified dispersion relations are found very easily by manipulating Eqs. (186|) . In momentum space 

A„(x) = j ' d A x A^e- 1 ^' (92) 

and choosing the Lorentz gauge, these equations reduce to 

k 2 A^ + 2 Pfi (p v A v ) = , P a = 2VBD a0 kf3, (93) 
k v A v = . (94) 

The vector p a is proportional to the momentum space version of the vector introduced in Eqs. ([55)1 and ([55")) . 
Multiplying the first relation in ([93]) by p^ , it follows that 

(k 2 + 2p 2 )(p»A v ) = 0. (95) 

Moreover p v ' A v is gauge invariant. In fact, in coordinate space is proportional to D a @ f a p . So this component is 
physical and has a dispersion relation given by 

k 2 + 2p 2 = 0. (96) 

If (k 2 + 2p 2 ) is not zero in Eg . ([9"5|) . it follows that p v A v = . In four dimensions, this condition plus the Lorentz 
gauge leaves two degrees of freedom. But the Lorentz gauge permits a further gauge transformation with parameter 
A such that 

d 2 X = 0, (97) 

which leaves only one degree of freedom as it should be. Moreover from the first equation in (|93p . this remaining 
degree of freedom satisfies 

k 2 A^ = 0, (98) 

so its dispersion relation is 

k 2 = 0. (99) 

The general solution of ([95|) is 

p"A v = -Ai(fc)^(fc 2 + 2p 2 ). (100) 

Putting it back into ([9"3"|) , we get 

A^ = a^k)5(k 2 ) + 2^^S{k 2 + 2p 2 )p^ a M (fc)fc" = 0, a M (fc)j>" = 0. (101) 
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From (jlOip . we obtain the electromagnetic tensor 

Uu = (k^a v (k) - k v a^k))5{k 2 ) + 2^±6(k 2 + 2p 2 )(k^ - k vPfl ). (102) 

It represents a plane wave with dispersion relation k 2 = (the magnetic and electric fields are the standard ones, 
perpendicular to k), plus a plane wave propagating in the direction k with dispersion relation k 2 + 2p 2 = 0. 
The fields for the second wave are 

Ej = fa = 2(k oPj - k jPo )^l, Bj = 2e jk i(k m )^l. (103) 

Notice that E is perpendicular to B, B is perpendicular to k, but E is not necessarily orthogonal to k. Moreover, 
this wave exists only if A; 2 ^ 0. 

B. 3+1 Decomposition 

An alternative way to determine the polarization of the plane wave solutions is by splitting Eq. (|86p into the 
components 0, i, which yields 

(- V 2 + 2bD' 2 ) A - d d l Ai = SBDoDjAj , (104) 
(d 2 A k - d k d A + dkdiAi) = &BD k (DjAj + D A ) . (105) 



We choose the Coulomb gauge 



The final equations are 



d r A r = 0. (106) 



(-V 2 + SBDl) A Q = -8BD DjAj, (107) 

(d 2 A k - d k doA ) = 8BD k {D.jA 3 + D A ) . (108) 

As in the usual electrodynamics A can be expressed as an instantaneous function of Aj 

Ao = - 8B (-v4*BD>) DoD ^- (109) 



The final vector equation leads to 



OB 

d 2 A k = (-V 2 8 kt + d^) { _ v2 + 8]3D 2f iDjAj, (110) 

where it is a direct matter to verify the Coulomb gauge condition. 
The electromagnetic fields are given by 

A = {A% (111) 
E = {Ei}, Ei = fa = d Ai - dtAo, E = -9 A-VA°, (112) 

B = {Bi}, B k = -^e klm f lm = e klm d l A m , e 123 = +1, B = V x A, (113) 

which reduces to 

E r = -iujA r - ik r A , B k = -ie k i m kiA rn (114) 
in momentum space with the notation k = {ki} in the three dimensional subspace. 
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C. Modified Maxwell equations 

They are 

V-E + 86(e- V)( B-b-E-e) =47177 , (115) 

Vx (B + 8Bb( B • b — E • e)) — --^-[E + 8Be( B ■ b — E • e)l = —J , (116) 

cat c 

V-B = , (117) 
9B 

c dt 



1 <9B 

VxE+-— = . (118) 



In the notation of electrodynamics in a medium we can introduce 

D = E+8£e( B b E e), H = B + 8£b(B b E e), (119) 



such that 



In components 



Particular cases are 



1 <9D 4-7T 

VD = Airp, VxH — = — J. (120) 

c at c 

1 <9B 

V-B = 0, VxE+- — =0. (121) 
c at 



A = (S i:j - 8Be l e )E + 8Be i b j B j , (122) 
Hi = (5 i:j +8Bb i b j )B j -SBhejEj. (123) 



bj = -> Di = (S i:j - 8Be i e j )E j , H i =B i , (124) 
e* = -> A = A, A - (<% + 8Bb i b j )B j . (125) 

The first case induces a modification in Coulomb law 

B = 0, V-(E-8Be(e-E)) =4np, E = -V$, (126) 

(v 2 - 8B M,) 2 ) $ = -4np. (127) 

VI. SUMMARY OF THE DISPERSION RELATIONS AND ELECTROMAGNETIC FIELDS 

Here we use the previous decomposition of the background field in terms of three dimensional components 
together with the 3 + 1 formulation described in subsection (IV-B), to search for plane wave solutions of the Maxwell 
equations. We assume the space-time dependence of any field to be proportional to e 4 ( k - x - wt ) } where k is the momen- 
tum of the wave and we choose the potential A in the Coulomb gauge. The notation is 

b = {6 m }, e={e m }. (128) 

The expressions of Eqs. (155)) , (|59")l in momentum space are 

{D k } = D = -j[we-kxb], (129) 
D = -ie-k. (130) 

There are two main cases 

(i) DjAj — 0: In this case the dispersion relation is ui = |k|. The fields are 

B = 7{(e ■ k)k - e + bxk} , (131) 
E = 7 {(b-k)k-b-exk} . (132) 
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(ii) DjAj : In this case the dispersion relation and the fields are 



'256£ 2 (k • (b x e)) 2 - 4(1 - 8Be 2 ){8B[(k ■ e) 2 - (k x b) 2 ] - 1} - 16Bk- (b x e) 
w = |k|- — — ^ , (133) 



2(1 - 823e 2 ) 

B = 7{k x (k x b) — ujk x e} 
E = j{ojk x b - uj 2 e + k(e • k)} . 

In both cases 7 is an arbitrary constant. For small 23, the dispersion relation (|133[) reduces to 



= Ikl 



1 + B ( 8k • (e x b) + 4(k x b) 2 + 4(k x e) : 



(134) 
(135) 



(136) 



The above expression clearly shows that the model is stable in the small Lorentz invariance violation approximation, 
where the quantities Be 2 , Bb 2 , Z?|e||b| are very small compared to one. 

The anisotropic velocity of light arising from the above dispersion relation is 



V k w = c(k) = k (1 + 8B (b 2 + e 2 ) - 46 ^(k x b) 2 + (k x e) 2 

+8B (e x b) -8B (b • k) b-8B (e ■ k) e . (137) 

Next we consider the separate cases where e 7^ 0, b = 0, e = 0, b ^ 0, and e ^ 0, b ^ 0. The results are: 



We find the dispersion relation 



A. Case I: b = 0, e / 

1. CASE (1-1) e ■ A ^ 0, D A ^ 

, 1 -8Se 2 cos 2 0, o 



1 - 8Be 2 



(138) 



which agrees with the one obtained from (|55|) . Here 8 is the angle between k and e. Assuming the corrections terms 
to be very small we can rewrite (I138P as 



loi = |k| (1 + 4Z?e 2 sin 2 
The electromagnetic fields in momentum space are 

(e-k) 



E = a 



e — 



B = — [k x e] , 



(139) 



(140) 



where a is a dimensionless constant. The following orthogonality relations are satisfied 

E-B = 0, k-E^0, k-B = 0, (141) 

together with the properties 

E x B - [k-8S(e • k) e] = m, (142) 



The ratio between the amplitudes is 



m 

IBI 



m • B = 0, m • E = 0. 



sin 2 9 + (1 - 8Be 2 ) 2 cos 2 9 
{l-8Be 2 cos 2 9) {\ - 8Be 2 )' 



(143) 



(144) 
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2. CASE (1-2) e ■ A = 0, D ■ A = 
The corresponding dispersion relation is 



"2 

in accordance with (|99|) , with the electromagnetic fields 



u\ = k 2 , (145) 



E = (k x ej , B = (3 k x [k x ej , (146) 
where (3 is an arbitrary dimcnsionlcss constant. Also we have 

|E| = |B| =p\e\sm6. (147) 

The orthogonality relations are 

k-E=0, k-B=0, E B = 0, (148) 

similarly to the standard case. Also 

IFI 2 

[ -X = 1. (149) 
|B| 2 

B. CASE II: e = 0, b / 0, D = 0, 
1. CASE (11-1 j b • A ^ D • A = 

We find the dispersion relation 

uj\ = k 2 . (150) 

and we have 

E = 7 kx(kxb) , B = - 7 (kxb,) (151) 



with 7 an arbitrary dimensionless constant. Also 

, -1, k-E = k-B = 0, E-B = 0. (152) 



IEI 2 



|B| 

2. CASE (11-2) b-A = ^-D-A/0 

Here the dispersion relation is 

ul = k 2 (1 + 8Bb 2 sin 2 6) , (153) 

with 9 being the angle between k and b. The above result agrees with that obtained from (155]) . The electromagnetic 
fields are determined, up to the arbitrary constant 5, by 

E = — (k x b) , (154) 

B = - kx(kxb) (155) 

and they satisfy 

k-E = 0, k-B = 0, E-B = 0, (156) 

|E|2 - W * = (l + 8£b 2 sin 2 tf). (157) 



IBI 2 Ikl 2 



Notice that half of the cases previously considered violate rotational invariance, with the exception of those having 
unmodified dispersion relations together with standard properties of the propagating electromagnetic fields. 
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C. CASE III: b 2 = e 2 . e b = 



We consider this case corresponding to a plane wave vacuum because the final theory remains invariant under the 
three-parameter subgroup HOM(2). 
We choose a reference frame such that 



n=(l, 0, 0), b=(0, 0, b), e=(0, sb, 0), s = ±1, 
k- A = 
DjAj = ib ((fci - us) A 2 - k 2 Ai) 



(158) 
(159) 
(160) 



1. Case D k A k = 0. 

The dispersion relations are the usual ones. From the master formulae (116-117), the fields are 

B = 7{(s(b x n) • k)k -s(bxn)lbxk) , 
E = 7 {(b -k)k- b + sk x (b x n)} . 



(161) 
(162) 



The exact modified dispersion relations are 

1 



0J 5 



(1 - 8Bb 2 ) 

In the small violation approximation Eq. (|163p reduces to 

w 5 = Ikl [ 1 + 4B6 2 



2. Case D k A k 0. 



(k 2 - 86b 2 (k 2 - k 2 )) + 8Bb 2 sh 



ki 



From the formulae (118-119), the fields are 

B = 7{k x (k x b) — W5sk x (b x n)} 
E = 7{w5 (kxb)- w|s(b x n) + sk((b x n) • k)} 



(163) 

(164) 

(165) 
(166) 



VII. THE MODEL AS A SECTOR OF THE SME 



In order to impose some bounds upon the parameters of the model it is convenient to recast the quadratic sector 
of the action (|38p in the language of the Standard Model Extension and to make use of the numerous experimental 
constraints derived from it ( see for example Ref. [HI) an d summarized in Table II: Photon-sector summary, of Ref. 
[29[ | . To begin with, let us recall the dimension of the fields and parameters involved in the model 



[A] = m, 



[b]=m 2 , [B] 



1 



Be 2 



0. 



Next we make the identification 



where the tensor (kp) K with 19 independent components, has all the symmetries of the Riemann tensor and a 
vanishing double trace. In terms of the matrix elements characterizing the vacuum expectation values of the 
electromagnetic tensor we obtain 



(167) 



(168) 



4B 



1 



(169) 
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where 



D z =D al3 D ap = 2{h z -e z ). 
Next we have to identify the appropriate combinations of the components of {kp 



step in that direction constitutes the definitions of the components (kde) J i {khbY an d (kdbY 
which are written in terms of our parameters e, b and B as follows 

(K DE y k = -2 (k F ) 0l0k = -l2Be jek - BD 2 S 



which are bounded. A first 

kj 



'jk, 



(kmeY 11 = ^ejp q ekrs {k F ) pqrs = \2Bbjbk - BD 2 S jk , 



(kdb) 



jk 



(K HE ) kj = e kpq (k F )°^ = l2B( ej b k - i (e • b) S jk ), tr 



jk 
i DB 



= 0. 



The final combinations in terms of which the bounds are presented turn out to be 



(*e-) jk 



^ (kde + KHB) Jk = 66 



bjb k - eje k - ^ (b 2 - e 2 ) S jk 



< 10 



-32 



(k de - K, HB y k - ^Sj k tr(K DE ) 



1 



-6B 



{b 3 b k + e 3 e k )--8, jk {h 2 +e 2 ) 



(k db + K HE y k = 6B(ejb k - e k bj) < 1CT 



< 10 



-16 



(KDB — K-HE 



6B 



(ejb k + ekbj) - - (e ■ b) S jk 



< 10 



-32 



Rtr = -tr(k DE y = -2B(e z + < 10 



-15 



(170) 

\ fi 

(khe) 

(171) 
(172) 

(173) 
(174) 

(175) 
(176) 
(177) 
(178) 



The bounds for n tr have recently been improved from 10 7 [2!|, to 10 11 [3(|, and finally to 10 15 [3l|. In quoting 
our results arising from the bounds upon (K e+ y , (R a ^y k and Rtr we have chosen our reference frame in such a way 
that only the light-cone gets modified by the LIV parameters, while the propagation of the fermioms is the standard 
one 

All matrices from (|174[) to (|177p are traceless, with (/to+) being antisymmetric (3 independent components) while 
the three remaining ones are symmetric (5 independent components each). The above combinations (|174|) to (|178p 



constitute a convenient alternative way to display the 19 independent components of the tensor (kp 



The 



bounds are obtained in terms of such combinations referred to the Standard Inertial Reference Frame defined in Ref. 
[2~ij and are written at the end of each of the above equations. 
Notice that the bilinears which are odd under the duality transformations 



b —>■ — e, 



(179) 



are much more suppressed than those which are even. In this way, even though our model is not duality invariant, this 
transformation seems to explain the above mentioned hierarchy in the LIV parameters exhibited in Eqs. (|174p ~ (|178p . 
In order to obtain some specific consequences of the bounds (|174[) to (|178[) it is convenient to express them in the 
coordinate system defined by (|4"Tj) . Also we introduce the notation 



x = VOBexlO 16 , y = VOBbxlO 16 , 



(180) 



We consider only absolute values of the related quantities and we focus upon the stringent constraints in Eqs. (I174p 
and (|177|) . The non-trivial contributions are 

(181) 
(182) 
(183) 
(184) 



K+) 




{x 2 -y 2 )\ <3, 


/- n22 

K+) 




(1 -3sin 2 V)x 2 -y 2 \ < 3, 


1- n33 




2y 2 + (1 -3 cos 2 '!/') x 2 < 3, 


/- n23 




x 2 sin 2-0 1 < 2, 
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FIG. 1: Boundaries of the allowed region obtained from the constraints in the parameters k{, + and H J k _. The allowed region 
is on the inside the dashed, dot-dashed and dot-dot-dashed lines. An excellent approximation for it is the sector of the circle 
shown in solid line [33l ]. 

together with 

(k _) U = (k _) 22 = \xycosip\ < ~, (185) 

(k q _) 33 = \xycoaip\ < -, (186) 
(k _) 23 = \xysin4>\ < 1. (187) 

The allowed region in the (x — y) plane is shown in Fig.l, where we plot the boundaries of the corresponding 
inequalities. In expressions (|184[) , (|185[) - (|187p we consider the lower bound corresponding to the maximum value of 
the trigonometric function on the LHS of the inequality. This leads to 

x < y/2, \xy\ < | < 1 < |. (188) 

The boundary xy = 3/4 is plotted in dashed line. The most stringent bound from Eq. (|183|) results when ip — n/2 
and corresponds to 



j/ < y — 2 — ■ (189) 

The boundary is shown in dot-dot-dashed line. The most stringent bound from (|182[) arises again from ip = n/2 and 
corresponds to 

y < V3 - 2x 2 , (190) 

which boundary is plotted in dot-dashed line. The expression (|18ip does not provide additional bounds and is plotted 
for completeness. The corresponding boundaries are y± — y/x 2 ± 3 shown in dotted lines in Fig. 1. An upper bound 
including all previous ones is given by the interior of the circle 



ir=\\\~x\ (191) 



shown in solid line in Fig. 1 and which translates into the restriction 



B(e 2 + b 2 ) = -ij(e 2 + b 2 ) < 2.5 x 1(T 33 . (192) 
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This bound incorporates all the constraints established in Eqs. (|174[) - (|178p . All the above relations are valid in the 
Standard Inertial Reference Frame centered in the Sun (24[. 

VIII. SUMMARY AND FINAL COMMENTS 

We study a novel way of implementing a model with spontaneously broken Lorentz symmetry by introducing a 
constant vacuum expectation value (VEV) of the field strength (F^) = C^ v instead of imposing a VEV of the 
electromagnetic potential (A^) as frequently investigated in the literature. In this way, our model preserves gauge 
invariance from the very beginning. We start from the effective Lagrangian (|13p describing non-linear electrodynamics 
and containing a potential V e f f{F^u), which is argued to arise after integrating massivegauge bosons and fermions in 
an underlying conventional theory. This approach is an extension of the model in Refs. [7|,|8j|, which includes all quartic 
fermion interactions allowed by the Dirac algebra. For simplicity we have considered only one fermion species. The 
subsequent integration of them produces the effective gauge invariant potential which is bounded in the high-intensity 
field limit and has an stable minimum. In order to investigate the dynamical consequences of the electrodynamics 
constructed around such minimum we have chosen the standard Ginzburg-Landau parametrization for V e f f{F llv ). We 
have verified that the extremum conditions in the potential corresponds to the extremum conditions of the energy 
arising from (|13p . under the requirement that the involved fields are constant at the extreme points. Next we expand 
the fields around the minimum and arrive to the spontaneously broken action ()25|1 from where we can further discuss 
the resulting non-linear electrodynamics, which differs from the standard one described in Ref. [26], for example. The 
action (|25|) contains the field strength a M „ (E , B), together with the auxiliary field X M which equations of motion 
demand that the two-form a is such that da = 0, thus introducing the electromagnetic potentials. The constitutive 
relations expressing the electromagnetic excitations D and H , are identified via the equations of motions (f2l)|) , (f3~lj) 
and provide a direct connection with the field strengths E and B. After eliminating the auxiliary field and writing 
the most general solution to the condition da = 0, we arrive at the action (|38|) . which is the starting point of our 
subsequent analysis. Some redefinitions allow the symmetry breaking parameters to be characterized by the constant 
antisymmetrical tensor (proportional to the VEV C^), plus an additional constant 6. 

The search for the subgroups under which the model remains invariant after the symmetry breaking includes the 
generator of dilatations plus the standard Lorentz generators. In this restricted subalgebra the former commutes with 
the remaining ones, which allows its realization as a multiple of the identity. There is only one case, given in the 
subsubsection IV.E-2, which incorporates this generator, resulting in the breaking of the Lorentz group to the three 
generators subgroup HOM(2). All the other cases break to a subgroup isomorphic to T(2), with two generators. 

The modified photon dispersion relations together with the corresponding plane wave polarizations are classified 
according to the value of Di a k a A l , where k a is the plane wave four momentum and A 1 is the corresponding vector 
potential in the Coulomb gauge. The amplitudes of the propagating fields are written in terms of those parameterizing 
D^ u . The case Di a k a A l = leads to dispersion relations u> = |k| with the triad E, B, k having standard orthogonality 
properties. The situation Di a k a A 1 ^ produces the following properties for the vectors involved: E is perpendicular 
to B, B is perpendicular to k, but E is not necessarily orthogonal to k. The corresponding dispersion relations are 
of the form 

uj ± = |k| F ± , (193) 

where F± is independent of |k|, being only function of the angles between k and e , k and b, k and e x b. In this 
way our model predicts anisotropy in the speed of light. Also, as shown in Eq. (|136[) the model is stable in the small 
Lorentz invariance violation approximation where the quantities 6e 2 ,66 2 ,6|e||b| are very small compared to one. 
This is validated by the bound 6(e 2 + b 2 ) < 2.5 x 10^ 33 derived from Fig.l. 

In order to make a more quantitative statement about this anisotropy let us consider the different possibilities 
arising from the velocity (|137p . To the leading order 6|e| 2 , 6|b| 2 , 6|e||b| the magnitude of such velocity is 

c(k) - 1 + 86 (e 2 + 6 2 ) - 46 (b ■ k) 2 + (e ■ k) 2 - 2k- (e x b)^ , (194) 

which we take as the one-way light speed in the direction k, predicted by our model. The corresponding two-way 
light speed is then 

c TW (k) = i (c(k)+c(-k)) = l + 86(e 2 + fo 2 ) - 46^(b-k) 2 + (e-k) 2 ^) . (195) 

A standard measure for the anisotropy of the speed of light is the comparison of the two-way velocities in perpendicular 
directions. We single out the vector (e x b) = n in our reference frame, which plays an analogous role to the relative 
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velocity among the ether frame and the laboratory frame in the Robertson-Mansouri-Sexl (RMS) type of analysis 
[1H , [35l ] . The vectors k and n form a plane in which we define the vector q perpendicular to k given by 



kx ( k x n 



k • n ] k — n. 



Next we calculate the two-way light speed along this perpendicular direction obtaining 



ek 



crw(q) = 1 + 8B (e 2 + b 2 ) - 4B (k • nj \b 2 (b • k 

An appropriate definition for the anisotropy of the speed of light in this model is 

Ac 



c 

Ac 



CTw(k) - CTw(q) 
2 



46 1 



kn 



b 2 b ■ k 



„2 / 2 



ek 



From the last expression we obtain the bound 

A 



— < 4Bsin' 0(b 2 + e 2 ) < AB(b 2 + e 2 ) < 10" 

c 



(196) 
(197) 

(198) 
(199) 

(200) 



according to Eq. (|192[) . where 9 is the angle between the vector n and k. The above anisotropy measures the difference in 
the two-way speed of light propagating in perpendicular directions in a given reference frame. The standard Michelson- 
Morley type of analysis can be made by measuring first the time difference along two perpendicular trajectories 
Li (propagating with Crw(k)), and L 2 (propagating with ctw(<1)) 



Li 



CTw(q) Crw(k)' 

Subsequently a similar measurement is made by rotating the apparatus by 90 degrees 

/ L2 L\ 



CTiy(-k) c TW (q)' 



Finally, the difference in optical paths 



At= \t-t'\ =4B(L 2 + L 1 )sm 2 e [b 2 (b • k 



(201) 



(202) 



(203) 



measures the change in the interference pattern of the Michelson-Morley experiment, according to our model. 

It is important to emphasize that the bound (|200|) is not related to the anisotropy in the speed of light induced by 
the passage from a preferred frame (where propagation is isotropic) to another frame moving with relative speed v 
and subjected to alternative methods of time synchronization. As such, it cannot be directly compared with the RMS 
ty pe b ounds appearing in the literature. The difference in two perpendicular two-way speed of light in the RMS case 



is 3£ 



5c 



( v 

\c T w{4> + tt/2) - Ctw{4>)\ = ~ 

V c 



(P-- + 5) cos 20 



(204) 



Here <fi is the angle between the photon direction and the relative speed v. Recent bounds upon the Mansouri-Sexl 
parameter {P — \ + 5) [36j are of the order of 10~ 10 , implying a RMS bound 



- < io- 16 , 

c 



(205) 



with the relative speed between the earth and the CMB reference frame, where the cosmic background radiation looks 
isotropic, having the value v = |v| ~ 300 km/s. 

Perhaps a more significant comparison of our results can be made with the work of Ref. [37] , where the Euler 
model of nonlinear electrodynamics [38| is used to predict the following isotropic bound of the one-way speed of light 



5c 



< 1.2 x 10 



-23 



(206) 
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In our case this would correspond to averaging over all angles in (|194[) , leading to 

6c 



— < SB(e 2 + b 2 ) ~ 2 x 1(T 32 . (207) 

c 



Let us emphasize that in all our estimations we are assuming that the involved references frames (CMB, Sun based, 
earth based) are concordant, in such a way that first order quantities in the Lorentz invariance violation parameters 
remain of the same order in all of them. 

Assuming that our background fields e and b might represent some galactic or intergalactic fields in the actual era, 
we obtain a very reasonable bound for the magnetic intergalactic field by assuming that the constant appearing in 
the action (f38|) corresponds to an energy density p 

p=\[(l- D 2 B] D 2 ~ \ (b 2 - e 2 ) , (208) 

which we can associate with the cosmological constant, since it would represent a global property of the universe. The 
fact that this constant is positive favors (b 2 — e 2 ) > 0, so that one can perform a passive Lorentz transformation to a 
reference frame where e = 0. Supposing further that this frame, which describes the intergalactic fields, is concordant 
with the Standard Inertial Reference Frame and taking the upper limit [39j 

|pa| < 1(T 48 {GeV)\ (209) 

we obtain the bound 

|b| < 5 x 10~ 5 Gauss, (210) 

which is consistent with observations of intergalactic magnetic fields. Let us recall that 1 Gauss — 1.95 x 
10~ 20 {GeV) 2 . It is amusing to observe that this bound on |b| translates into the following condition for the 
mass scale Mb introduced in Eq. (|192j) 

M B > 1.4 x 10~ 4 GeV. (211) 

Recalling that the lightest charged particle is the electron with m e = 5.1 x 10~ 4 GeV , we notice that the above bound 
is consistent with the QED Euler-Heiscnbcrg type of non-linear corrections to the photon interaction, which scale as 
1/to 4 , with m being the mass of the charged particle contributing to the loop in the effective action [4fj |. 



The Appendix 



In this Appendix we extend the models in Refs 0, [1| to include all quartic gauge invariant fermionic interactions 
allowed by the Dirac algebra and look for their contribution to the effective gauge invariant electromagnetic action 
by integrating the fermion field. We heavily rely on the discussion and results of Ref. 25], which presents a complete 
and detailed description of the standard electromagnetic calculation. We do not attempt to cite the original papers, 
which references can be found in Ref. [25| . 

After integrating the massive gauge bosons in an underlying conventional theory Q, let us consider the effective 
Lagrangian 



L e t f = *(i'f (d^ + ieA^ 



2A 2 

M,a 



[(#M a #) (*M°#)] 



(212) 



The index M : S,V,T, PV, PS labels the tensorial objects that constitute the Dirac basis: Scalar, Vector, Tensor, 
Pseudovector and Pseudoscalar, respectively. The generic index a labels the covariant (contravariant) components of 
each tensor class. More details are given in Table I. We follow the conventions of Ref.[22j], with appropriate factors 
chosen in such a way that each current (\&M a \D r ) is real. 

Our goal is to estimate the contributions of the additional quartic fermionic couplings to the standard effective 
electromagnetic action W E X M . The latter is obtained by integrating the quadratic contribution of the fermions in 
Eq. (|212|) and it is given by 

<0 + |0_>fi) = 



[D#] [£>*] exp 



det(* {iY {dft + ieA^) - m) = exp iW^ M = exp 



d 4 x L 



(i) 

EM 



(213) 
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To this end we introduce the gauge invariant auxiliary fields (Cm)o,, via the couplings 

[Pa 
I 2 



L eff = § (i 7 " (fy + ieA^) - m) * + ^ [M(C M ) a (C M )° - 9Af(C M )a (*M°*) 



(214) 



M,a 



The equations of motion for the auxiliary fields are (we do not use the Einstein convention for repeated indices unless 
explicitly stated) 



which reproduce (1212[) when 



(C M )a = — (*M a *) 

Pm 



A 2 pm ' 



(215) 



(216) 



Next we integrate over the auxiliary fields following the standard steps which begin with the introduction of the 
currents (jm) together with the replacements 



(Cm) a 



iS{j M ) a 

The total vacuum transition amplitude can then be written as 



(0+|0- 



[DV] [£>*] [D(C M )a] exp 



/ / d A x L e ff 



exp 



I d ' X UM) a (C M )a 



(217) 



(218) 



with 



L eff = * | vf (d, + ieA,) -m-J2 mM a — U + £ \ p Jf(C M ) a {C M T 

M,a l °K3M) / M ± 



(219) 



Now we can perform the Gaussian integral over the fields (Cm) a, which includes the linear contribution of the currents, 
to obtain 



(0+|0_) A = J [£>*][£>*] exp i J d i x^\i^{d^ + ieA^)-m-Y^MM a ^^ a 

f d i xY / —(jM) a (jM)a , 
J PM 



x exp 



(220) 



where the limit (jM)a — > is implicit. Notice that in order for the functional integral of the fields (Cm) a to be well 
defined in the usual Euclidean analytic continuation we must require that p a > 0. 
Using the standard formula (See for example Eq. (1.53) in Ref. (2d| ) 



F 

















F(X) 


j exp 






exp 




exp 


\-- [ —A— 
2 J 5X 6X_ 



(221) 



X=A 3 

where A(a; — x') is proportional to 5(x — x') in each case, we can rewrite the above equation as 

5 S 



(Q+\Q-) A = exp 



_ z / d 4 xY — 
J ^ ^p> 

J M,a 1 



2p M S(X M ) a S(X M )a 



i j d A x * (i<f (d^ + ieAJ ~m)-J2 qAtM a (X M ) a * 

\ M,a J 



(X M )=0 



exp 



2 PM 5(X M ) a 5(X M )a 



det (*y (d„ + ieA„) - m) - ^ q M M a (X M )a 



M,a 



Xm=0 

(222) 
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where we have already taken the limit (ja/ ) a = 0. The next step is to rewrite 

(i^ (fy + ieAJ ~m)~J2 QMM a {X M ) a = (i-f (0„ + ieA^) - m) I 1 - G A ^ q M M a {X M ) a , 

M,a V M,a ) 

where G A is the fermion Green function in the presence of the external electromagnetic field A^, satisfying 

{iY (d^+ieA^-m) G A = 1. 

Also we recall the expression 

detP = exp [Tr In P] , 

where the trace is understood both in coordinate as well as in internal spaces. In this way we have 



(223) 

(224) 
(225) 



<0+|0_) 



exp 



iW 



(i) 

EM 



exp 



1 



M.< 



2p a 5(X M ) a 6(X M )a 



x exp 



Tr In [l-G A ^2q M M a (XM)a 



Ma 



(226) 



To make an estimation analogous to the two-loop corrections W^l 1 (A) to the effective action W^(A), which is 
already calculated in Ref. (25| . we expand the third exponent up to terms quadratic in (Xm)o' Calling the operator 

G A lMM a (X M ) a =G = J2 Sm (227) 



M.a 



M 



we are left with 



(0+|0_) A =exp 



iW 



(i) 

EM 



exp 



1 



M.a 



2p M 5{X M Y S(X M )a 



exp [iTr (i0)] exp 



-Tr (iQ' 



(228) 



x M =o 



With the purpose of having a preliminary order of magnitude estimation of the corrections we next consider separately 
each of the terms Qm in the summation appearing in Eq. (|227p . In this way we are neglecting the contribution of 
the interference terms among different currents j'm, even though their contribution could also be calculated along the 
lines presented here. 

In order to proceed we make use of another well-known identity 

~Tr In (1 - BAy 1 
- l ~ J X B{1- AB)- 1 X 
+i J X(l - BA^F 
+ ^ J FA(1 - BAy 1 F , (229) 





■ i r 6 _5J 








exp 


2JSX SX_ 


exp 


% - J XBX + i J FX 


= exp 



x exp 
x exp 
x exp 



for a fixed M. Next we identify the corresponding operators 



(A M y c (x-x') 



1 

Pm 



5°6 4 (x - x'), F Ma = iq M G A (x, x')M a , 



(B M ) c (x,x') = i q M G A {x,x')M G A {x',x)M c , 
where 8 b c denotes the identity in the corresponding case. For example 



M a - y, 8\ - 61 ; M a - a*", S b c -> \ (s^ p - S^ a ) 



(230) 
(231) 

(232) 
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M 


M c 


9m = Y,M MCM *) 


03M =J2r tr ( a 3M c a 3 M c ) 


S 


I&X& 


4 


4 


V 


7" 


16 





T 


a"" 


48 


-16 


PV 




-16 





PS 


75 


-4 


-4 



TABLE I: General notation and traces for the different contributions to Im in Eg. (|238p . 



In the sequel we also restrict ourselves to the case of a constant (or slowly varying) external electromagnetic field. 
In this case G A (x,x') is independent of position so that the transformation properties under the Lorenz group imply 
that F Ma = 0. (See for example Eq. (7.8) of Ref. [H|]). This further simplifies Eq. l|2"2"gl) to 



(0 + |0_)i=exp 



iW 



(i) 

EM 



exp 



-Tr\n{l- B M A M ) 1 



exp 



iW, 



(i) 



EM 



exp 



-Tr (B M A 



M , 



This defines the sought correction 



exp 



iW {2) 



exp 



+ -Tr(B M A M ) 



M fixed 



where 



Tr(B M A M )=i ^Y,J d4x tr{G A (x,x)M c G A {x,x)M c ) 



(233) 



(234) 



(235) 



where tr is the trace in the Dirac-matrices space. In other words we have 



u .,., J d 4 xtr(G A (x,x)M c G A (x,x)M c ) = ^ Id 



r (2) _ 2r M 

M ~ A2 



4 xl 



M- 



(236) 



In order to estimate each contribution we further consider the case where the external field is a constant magnetic 
B field in the z-direction. We avoid the constant electric field case because of the inherent instability due to pair 
creation in the strong field regime. For the case under consideration we have [ill ] 



G A (x,x) 



167T- 



f -~ exp [— ism 2 ] — e lCT32: , s — eBs, <r 3 = ij 1 "/ 2 

Jo s z smz 



(237) 



where m is the mass of the integrated fermion. In this way, the general expression for the quantity l^j defined in Eq. 
(f236]) is 



Im = 



ds ds' 



exp [—ism 2 ] exp [— ism 2 ] 



smz smz' 



(l6^ 2 ) J 

x ^2 [(cos z cos z'tr (M°M C ) - sin z sin z'tr (a 3 M c a 3 M c )) + sin(z - z')tr(a 3 M c M c )] 



(238) 



The last term in the second line of Eq. (|238p vanishes because ^ c M C M C is always a multiple of the identity and 
tr((T3) = 0. Let us introduce the notation 



9m = tr M c M^j , 6 3M = tr ^MVsM^ , 



which values are given in Table [U Also it is convenient to define 



^ r HI ^COSZ dS r . J-, dS r jn 

/ — exp \ —ism , 1 2 = / — tt exp — ism z = eB / — exp \ —ism 

Jo smz Jo s Jo s 



h = 



In terms of the above quantities we have 



I'M 



(l^s) [9mI1 ~ 9aMl - 



(239) 



(240) 



(241) 
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The second contribution, with infinite coefficient, behaves as B 2 so that we consider it as part of the renormalization 
procedure that we demand of the low energy limit of the effective potential 



hm V eff (B)=-pB 2 



p>0. 

In this way we consider the full expression for the contribution due to the additional currents to be 



(242) 



L 



(2) 



J2 L 



(2) 
M i 



(2) 
J M 



M 



A 2 



The effective Lagrangians in (|243p lead to the effective potentials 



y(2) _ 
V M — 



r (2) _ _ / m 

~ V167T 2 



(243) 



(244) 



The main point to be made is that each contribution to the effective potential arising from the subspace M contributes 
with an specific sign, as can be seen from Table 1. In these way, as it will be argued in the sequel, it is possible to 
have a net positive contribution for the total effective potential in the limit B/Bq — > oo, where Bq = m 2 /e is the 

(2) 

critical value of the magnetic field. Let us emphasize that Vf^ is still not well defined because the integral I\ diverges. 
This means that we should regularize it and further implement a renormalization prescription for the full effective 
Lagrangian 



L 



eff 



aB 2 



L%(B) + L^(B) 



(245) 



(2) 

Here L EM {B) refers to the two- loop standard electromagnetic correction. Our renormalization conditions will be 

fim L eff =pB 2 , p>0, (246) 

in such a way that the effective potential is a decreasing function in the vicinity of B = 0. The final normalization to 
the value of £ e // = —b 2 /2, 6^0, will be imposed at the end of the procedure, after we expand the magnetic field 
around the expected stable minimum Bmiu of V e f / , such that 



B = B 



Mm 



+ b 



(247) 



and we are left with the physical component b. The renormalization procedure should be analogous to the one carried 

out in Ref . [25| for the contribution L^ M (A) and it is rather involved due to the complicated field dependence of the 
regularized version of the divergent integrals. This calculation is out of the scope of these preliminary estimations. 

Thus we will only consider a particular regularization of I\ to proceed with the comparison of the different high- 
intensity field contributions to the effective potential. We take 



h = r exp 



-ism 



1 



(248) 



which produces a zero contribution when B — > 0. This integral can be readily calculated from the first order 
electromagnetic correction of the effective action given by 25] 



(i) 

J EM 



1 

8^ 



[°°ds , ai 
/ —5- exp —ism 
Jo s 3 



z cos z 1 

— + rz- 

sin z A 



1 



The result is 



h 



g T 2 ulj EM 



eB 2 
3m 2 



(249) 



(250) 



which yields our final result 



V M — 



16^ 



TM 

A 2 



1M 



2 ulj EM 



eB 2 



dm 2 3to 2 



(251) 
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In order to produce the required results in the low (high)-intensity field approximations, we list here the behavior of 
the relevant terms obtained from Ref. [25j. For B->oowe have 



V, 



(i) 



EM 



v, 



(2) 



24tt 2 

„4 



g ?r 2 ulj EM 



+ 



EM 



eB 2 



1287T 4 



B In 



-B 2 In 



- 



171 J 

eB 
m 2 



dm 2 3m 2 



eBhx 



2eB 



The limit B — > produces 



V, 



(i) 



BT {1) 
8tt 2 em 



EM 

eB 2 



-B\ 



dm 2 3m 2 



3607r 2 m 4 
( eB 2 1 e 2 5 4 
1,3m 2 " + 18 m 4 



In this way, the final contributions of a given current j'm to the effective potential are 



S 



oo 



V, 



(2) 



M 



(2) 
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M 



16tt 2 
1 



16tt 2 



(^) 2 r M e M e 2 B 2 
eB 2 
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3m 2 



/2e5\ 

1 e 2 B 4 
18 m 4 



The total contribution from the currents jm is 



(2) 



(252) 
(253) 
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(255) 
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(257) 
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